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We calculate the static and dynamic single-particle correlation functions in one-dimensional (ID) 
trapped Bose gases and discuss experimental measurements that can directly probe such correlation 
^^j", functions. Using a quantized hydrodynamic theory for the low energy excitations, we calculate 

■ both the static and dynamic single-particle correlation functions for a ID Bose gas that is a phase- 

' fluctuating quasi-condensate. For the static (equal-time) correlation function, our approximations 

, and results are equivalent to those of Petrov, Shlyapnikov and Walraven. The Fourier transform 

of the static single-particle correlation function gives the momentum distribution, which can be 
measured using Doppler-sensitive Bragg scattering experiments on a highly elongated Bose gas. We 
show how a two-photon Raman out-coupling experiment can measure the characteristic features of 
the dynamic or time-dependent single-particle correlation function of a ID Bose quasi-condensate. 
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I. INTRODUCTION 



In low-dimensional uniform Bose gases, the single-particle correlations are fundamentally different from a three- 
| dimensional gas, dramatically affecting the nature of the Bose condensation phase transition |Q. Gorlitz et. al. Q have 
created one- and two-dimensional trapped Bose gases by making the level spacing of the confining potential (along one 
or two directions) larger than the mean-field interaction energy. Recently, several experimental and theoretical studies 
have highlighted the manner in which the unique correlations of a uniform low-dimensional Bose gas are modified in 
a trapped Bose gas. With the increasing interest in BEC matter wave interferometry and atom lasers, it is clearly 
y—i \ important to have a greater understanding of one-dimensional (ID) trapped Bose gases. 

^ ■ In uniform systems, it is well known that the behavior of Bose gases is quite different in one- or two-dimensions com- 
pared to three-dimensions. In 3D, below a critical temperature Tbec, the long-range correlations of the phase are man- 

' ifested in the asymptotic properties of the single-particle correlation function, where limir^oo {r)^){Qi) \ ~ no ^ 

O^l ■ indicates long-range order in the system. Above TbeCi this single-particle correlation function decays exponentially 
for large distances, characteristic of normal systems. In systems of reduced dimensionality, long wavelength fluctu- 
ations of the phase of the Bose order parameter lead to a new state intermediate between a Bose-condensed and 
non-Bose-condensed system. The system forms a quasi-condensate, in which the phase of the order parameter is only 
coherent over a finite length scale. The system may be thought of as Bose condensate sub-systems which are well 
defined locally but are not correlated with each other M . 

The presence of a quasi-condensate causes the equal-time single-particle correlation function to spatially decay with 

a characteristic power law, limi,.^^ ^^(r)^(0)^ ~ 1/ |r|", where the power law parameter a is characteristic of the 

system. In a uniform two-dimensional (2D) Bose gas, a Bose condensate can be shown to only occur at T = 0. At 
, low but finite temperatures, however, a 2D Bose gas can form a quasi-condensate. In a one-dimensional (ID) uniform 
Bose gas, true Bose condensation never occurs, but at T = a quasi-condensate is possible. The central point is that 

. ,— i \ the characteristic features of ID and 2D quasi-condensates are associated with long wavelength, low frequency (i.e., 

^ ■ long distances and times) phase fluctuations of the Bose order parameter. 

In a trapped Bose gas, the effect of the long wavelength phase fluctuations is complicated by the finite size of the 
system. While a quasi-condensate is still locally coherent, the size of the system will play a crucial role. Recent 
theoretical work on trapped Bose gases has focussed on the static correlation properties of trapped ID and 2D Bose 
gases as a measure of the long-range phase coherence of the Bose gas. Early work by Ho and Ma (3j considered 
primarily the zero-temperature (vacuum) phase fluctuations. They showed that the strong phase fluctuations were 
restricted by the finite size of the system and as a result the system could form a Bose condensate even at T = 0. Since 
the low lying collective modes have finite energy in a trapped Bose gas, at sufficiently low temperatures the thermal 
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fluctuations are frozen and the T = analysis still applies. The role of thermal fluctuations in determining the phase 
coherence of the trapped Bose gas has been recently discussed in considerable detail by Shlyapnikov and coworkers 
H Hi §|> 01- They demonstrated |j| that below a certain temperature T^, the long wavelength phase fluctuations in a 
trapped ID Bose gas arc effectively frozen. As a result, the quasi-condensate of the uniform gas becomes a true Bose 
condensate in the trapped gas for T < T$. This work was later extended to deal with the effect of phase fluctuations 
in a highly elongated (cigar-shaped) 3D Bose gas ||, which was shown to behave as a quasi-lD Bose gas. It was 
found that the long wavelength phase fluctuations above a certain temperature T% D would still destroy long-range 
order, leading to a quasi-condensate. 

The role of phase fluctuations in an elongated 3D Bose gas has been recently explored by Dettmer et. al. ||, where 
evidence was found for a phase-fluctuating quasi-condensate by observation of characteristic interference effects in the 
freely expanding gas. This kind of real space interferometry experiment is a measure of the equal-time single-particle 
correlations of a Bose gas. However, the theoretical analysis required to evaluate the effect of phase fluctuations is 
quite involved and somewhat indirect Q . The Doppler-sensitive Bragg scattering experiments recently performed in 
Orsay || give an alternative and perhaps a more direct probe of the single-particle correlations, albeit in momentum 
space. 

Such Bragg scattering techniques have been used to measure the excitations in a trapped 3D Bose gas |nj [Ll| [l2|, [l^, 
0. Out-coupling (or tunneling) experiments discussed in Refs. [l6| [lTj also show how one can obtain information 
about the dynamic single-particle correlation functions of a trapped Bose gas. There has been some related theoretical 
work on the absorption line shape of a uniform ID Bose gas p8| . This type of experiment is also similar to tunneling 
between quantum wires used to probe the single-particle correlations in ID electron systems [[l9[ p0[ . 

We restrict ourselves to the trapped ID mean-field Bose gas, as defined in Refs. @, fj], ^3J. In Section 0, we 
summarize the calculation of low energy collective modes responsible for the long wavelength phase fluctuations using 



a quantized hydrodynamic theory |22|, |23|. In Section III, we use these results to calculate the static correlation 
function for a quasi-condensate, following the work of Shlyapnikov and coworkers |5|, ^, [7]]. The results for the 
static correlation function of a quasi-condensate are compared with that of a true Bose condensate in the ID trap at 
T = 0. We then calculate the momentum distribution of atoms both in quasi-condensates and true Bose condensates, 
since this is the key quantity measured in Doppler-sensitive Bragg spectroscopy. This discussion of results for the 



static correlation function sets the stage for our generalization to ID dynamic correlation functions in Section IV. In 
Section we show how this dynamic correlation function could be measured using a two-photon Raman scattering 
(out-coupling) experiment. 

II. CORRELATION FUNCTIONS 

The two-photon optical Bragg scattering cross-section can be shown to be proportional to jl2[ 

[S(q,w)-S(-q,-w)], (1) 
where the dynamic structure factor is defined by 

S(q,u) = J dte lut J dr J dr'e"^'^- 1 "') (n(r, i)n(r', 0)) , (2) 

with n(r) = t/>T (r)^(r) being the local density operator. Quite generally, one can prove from detailed balancing that 
S(— q, — lu) = e~^ u S(ci, uo). Even for a trapped gas, we can express the quantum field operators in the momentum 
representation 

1 ik-] 

in which case (0) reduces to 



S(q,u) = — jj^j^ <lfc' ui (^ t (k,^(k + q,iW t (kW-q))- (4) 

As we review in an Appendix, within the impulse approximation (IA), one can reduce Eq. (|i|) to the simple looking 
result 



/dk 
-^—^n(k)5(uj - q 2 /2m - q • k/m), 



(5) 
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where n(k) = (?/'^(k;)'0(k)) is the thermally averaged equilibrium momentum distribution of atoms in a fully interacting 
trapped gas. This result also shows that the energy transfer in the I A is very large and hence the second term in 
Eq. ([!]) can be neglected (w > k B T). 

Thus one sees how Doppler-sensitive Bragg scattering gives information about the atomic momentum distribution 
n(k). In turn, this momentum distribution is the Fourier transform of a single-particle correlation function, which we 
now discuss for a ID Bose gas. We consider a Bose gas in the mean-field regime which is trapped in a highly elongated 
cigar-shaped trap. In this ID geometry, the confinement along the radial direction is so tight that the trap frequency 
along the radial direction u p is much greater than the mean-field interaction energy. At the low temperatures of 
interest, the dynamics of the atoms in the radial direction are essentially "frozen", with all the atoms occupying the 
ground-state of the harmonic trap in the radial direction 




-(xW)/^ 



(0) 



Here the extension of the wavefunction in the radial direction is given by the harmonic oscillator length l p = yjh/rruv p . 

In the mean-field regime, l p is much greater than the radius of the interatomic potential R e . The scattering of 
atoms in this ID system is thus still three-dimensional. According to Ref. p5[ , the effective coupling constant in this 
ID system is 

2h 2 a 

9id - 777 ttt, (7 ) 

m l p (l p — La) 

where a is the s-wave scattering length and C is a numerical constant on the order of 1. The term Ca in Eq. ([?]) is 
negligible for lg_p> R e - In this limit, the expression for g\o is the same as that obtained by averaging over the radial 
wavefunction (Q) |f| j25fl, 



(8) 



9id = g dxdy4> p {x,y) 

J —oo 

2h 2 a 
ml 2 



We use this expression in the rest of this paper. 

Since the trap frequency along the radial direction is much greater than the mean-field interaction energy, in the 
mean-field regime we can construct a ID Gross-Pitaevskii equation for the axial direction only. Within the Thomas- 
Fermi regime, the ID Gross-Pitaevskii equation leads to a density distribution in the axial direction given by 



where the chemical potential is 



no(z) = ~U l0 ~ \ muJ l z2 ^j 6 (mo - \rau? z z 2 j . (!)) 



'A^f- ( io ) 



,4-\/2 hw z l 

The Thomas-Fermi length of the atom cloud in the axial direction is then defined by 



W^- (ID 



The number of atoms is given by N and l z = \Jhjmuj z is the harmonic oscillator length in the axial direction. 

The theory of the correlation properties of a phase fluctuating Bose gas builds upon this zero temperature theory 
of the static condensate density. The central idea is that the time-dependent single-particle correlation function 

C^(ry;i) = (^(r,^(r')} (12) 

can be calculated using the T = low energy collective modes calculated as fluctuations about the static condensate 
density, even when we arc dealing with a quasi-condensate with large phase fluctuations at finite temperature. In cal- 
culating the single-particle correlation function in Eq. (ft2|), we write the Bose field operator as -0(r, t) — yjn{r, t)e l ^ r '^ 
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using the density (n) and phase (</>) variables. While the density fluctuations are small, the phase fluctuations can be 
large in ID and 2D Bose gases. Within this approximation, only the phase fluctuation dynamics are considered and 
then the single-particle correlation function is given by 

C^(r,r';t) ~ y/n (r)n (r') (e~i^-^) . (13) 



It is adequate to use the static condensate density no(r) at T = given in Eq. (|13J), since this is unaffected by phase 
fluctuations at low temperatures. 

To evaluate the single-particle correlation function in Eq. (^), we need the normal modes corresponding to low 
energy phase fluctuations. We use the quantized hydrodynamic theory developed by Wu and Griffin p3[ . The 
fluctuations in the density and phase about their equilibrium values are given by 

h(z,t) = n (z) + Sh(z,t) (14) 

and 

fa,t) = Mz)+8fc,t). (15) 

By starting with the Gross-Pitaevskii Hamiltonian, we may develop a quantized hydrodynamic Hamiltonian describing 
the low frequency excitations at T = 0. Working to lowest order in the fluctuations and making the Thomas- Fermi 
approximation, one obtains the effective hydrodynamic Hamiltonian |23j ] 

/oo 
dz [mno(z)5v z (z) + g\D$n ; (z)] , (16) 
-oo 

in which the fluctuations of the condensate phase and density are treated up to quadratic order. The velocity 
fluctuation operator 5v z is related to the phase fluctuation operator in the usual way, 

5v z (z,t) = -^54>(z,t). (17) 
m dz 

The quadratic Hamiltonian in Eq. ( |l6|) has the same form as the renormalized ID hydrodynamic Hamiltonian derived 
by Haldane (^6), which has been used in recent studies on ID Bose gases p7| . 

Since the hydrodynamic Hamiltonian in Eq. ( |ilf ) is quadratic in the density and phase fluctuation operators, one 
can diagonalize it using the canonical transformations |23| 



Sn(z, t)=J2 ^jWe^'^j + A*^*(z)e^ t a]] (18) 



and 



54>(z, i)=J2 [ B j^( z )e~ iUJjf aj + B j ^We^'aj] ■ (19) 

3 

This form of the canonical transformations ensures the Hermiticity of the density and phase operators. The operators 
ctj and a- annihilate and create excitations with energy hujj, and satisfy the commutation relations [ojjOj/] = 
[dj,o^,] = and [dj,c^,] = The mode expansion coefficients Aj and Bj are to be determined along with 

the mode eigenfunction ipj( z ) and the excitation energy hiOj psfl . The mode eigenfunctions are assumed to form a 
complete, orthonormal basis. From the equal-time commutation relation [Sh(z), S(j)(z')] = iS(z — z'), the coefficients 
Aj and Bj may be shown to satisfy the relation A*Bj = —i/2. 

Using the Hamiltonian in Eq. ( |l6| ) , the density and phase fluctuation operators satisfy the Heisenberg equations of 
motion 

d -^ = -§- z [n (z)6v z (z,t)] (20) 



and 

dSv z (z 7 t) d 
dt dz 



\gi D 6n(z,t)]. (21) 



5 



These two equations may be combined to give the decoupled equations 



d 2 5n(z,t) d 
dt 2 = di 



and 



d 2 S4>(z,t) 

dt 2 



d_ 

c)z 



^-n (z)-^-Sh(z,t) 
m az 



—no (z) ^-64>{z, t) 
m az 



(22) 



(23) 



Using Eq. p0| ) or Eq. (21) together with the relation A*Bj — —i/2, the expansion coefficients Aj and Bj are given by 

(24) 



2~g~ 



ID 



B 3 = 



9lD 



Using Eq. (E31), the equation for the eigenfunctions ipj is 



9id 9 
m dz 



d 

no(z)—il)j(z) 
az 



+ uj 2 i; J {z) = 0, 



(25) 



with the static density no (z) given by Eq. ([)]). Eq. ( p5| ) then reduces to a Legendre equation, with the eigenfrequencies 
given by ujj — uo z y/ j{j + l)/2 p8[ . The associated eigenfunctions ipj are given by 

iPj(z) = J 3 -±^Pj(z/Z ), 



(26) 



where Pj is the usual Legendre polynomial and Zq is defined in Eq. (|1 1|) . This expression for ipj is normalized to 
unity and the integer j runs from 1 to oo. Using these eigenfunctions, the density and phase fluctuation operators in 
Eq. |J) and Eq. © are 



and 



Sh(z,t) = ^2 

3=1 



S$(z,t) = J2 
J'=l 



hujj j + 1/2 



2.9 



ID 



Zn 



Pj(z/Z )e- iu '' t aj +H.c. 



(27) 



9iD 



2hiUi 



Z 



Pj (z/Z )e- iu, " t a j +H.c. 



(28) 



Identical results are obtained in Refs. |], ^] by directly solving the coupled Bogoliubov equations for the ID Bose 
gas. We note that the low energy modes of a highly e longated tr ap are very similar to those of the ID Bose gas 
considered here. In the former case, we have iOj — t^zyjU + 3)/4, while the wavefunctions are Jacobi polynomials 



instead of Legendre polynomials 



As noted in Ref. ||, for low temperatures, the density fluctuations described by (5h(z)8h(z')) are small in compar- 
ison with the phase fluctuations (6(p(z)S4>(z')) . We can understand this from the fact that for reasonable parameter 
values, we have g±o > fiWj for the low-lying normal modes, which in turn implies Aj <C Bj. This justifies ignoring 
density fluctuations in reducing the correlation function in Eq. (|lj) to Eq. (|l3). 

Within this mean-field approximation for the phase fluctuations, the dynamic single-particle correlation function 
in Eq. (Q3h is given by 



Cy,t,/,(r,r';i) = <j> p (x, y)4> P {x' ', y')y/n (z)n (z') 



(29) 



Here we have used the standard results for functional integration jlj involving a quadratic Hamiltonian in reducing 
(exp(— i[5(f)(z, t) — 6<t>(z' , 0)])) = exp(— t) — 5<p(z' , 0)] 2 )). For notational convenience, we define F(z,z';t) = 

([6<p(z, t) — 54>{z' ', 0)] 2 ). In Section |TJ, we evaluate and review the properties of this single-particle correlation function 
CL,t-^(r, r'; t) in the static limit (t — 0). Section IV considers the extension to non-zero times. 
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III. STATIC CORRELATION FUNCTION 



Up to now, only the static correlation function has been investigated in a trapped low-dimensional Bose gas. 
Dettmer et. al. 11 have provided impressive interferometric evidence for phase fluctuations in a tight cigar-shaped 
3D Bose gas. Doppler-sensitive Bragg scattering experiments provide an alternative probe of the static single-particle 
correlation function since they measure the momentum distribution of atoms in the quasi-condensate |j] . In such a 
Bragg scattering experiment on a quasi-lD highly elongated 3D system, a large momentum kick q is given to the 
atoms by a Bragg pulse directed alon g th e long axis of the trap. The signal is proportional to the dynamic structure 
factor in the high momentum regime [Q, the I A expression given in Eq. (||). 

The results of the present section are equivalent to those of Shlyapnikov and coworkers p, ra, fjj. However, in 
previous studies, the results were used to discuss interference patterns of expanding gases j^j. In this section, we work 
directly with the dynamic structure factor measured in Bragg scattering experiments. These results will be needed 



in our generalization in Section [V to time-dependent single-particle correlation functions in ID Bose gases. 

As in the case of a highly elongated 3D Bose gas, the dynamic structure factor is the experimental signal in a 
Doppler-sensitive Bragg scattering experiment on a ID Bose gas. For the two Bragg laser beams oriented so that 
the momentum kick given to the atoms is along the long axial direction, the dynamic structure factor in the impulse 
approximation given by Eq. (||) is 



dk 

Sia(<1 = <?Aw) = / -^—^n(k)6(u> - hq 2 z /2m- hq z k z /m). (30) 



This may be reduced to 



Sia{<1 = Qz2,uj) = j- ,^ f- — / dk x dk v n(k x , k y ,k z = y/ m/2frhj r {oj - U3 r )), (31) 

(27f) V 2r ^r J-oo J-oo 

where we have introduced the recoil energy huj r = h 2 q 2 /2m. This shows how the IA dynamic structure factor is 
broadened as a result of the Doppler effect, which gives a fairly direct measure of the width of the momentum 
distribution of the atoms. The standard approach to estimate the broadening of the dynamic structure factor is to 
fit the data to a Gaussian centred at the peak value u> = u> r , 



SiA{q, to) = SiA{q, to = w r ) exp 



(uj — uj r y 



2A 2 



(32) 



where Aia is the Doppler width of the dynamic structure factor, given by 

Aj A - - [S IA (d, u)/d*S IA (ci, u)] \ u=Ur . (33) 

For further discussion, see Ref. [ pi] . 

We now turn to evaluating the dynamic structure factor in Eq. ( |3l| ) in two different regimes: a Bose-condensed and 
a quasi-condensate gas. We first calculate the momentum distribution. Since both regimes occur at relatively low 
temperatures, we ignore thermally-excited density fluctuations. However, since these do not significantly affect the 
transition from a Bose condensate to a quasi-condensate as the temperature increases above [see Eq. ([I3|)1, we are 
not neglecting any key features that a Doppler-sensitive Bragg scattering experiment probes. A more sophisticated 
theory must be used to take into account the density fluctuations that become increasingly important at higher 
temperatures (see, e.g., Refs. |29|, pop). 

We first consider the momentum distribution for a true Bose condensate at T = 0. The static correlation function 
can be split into a static Bose condensate term and the density fluctuations about the condensed state 

(^(iW)) - $S(r)$o(r') + (V^(r)^(r')) . (34) 



The density fluctuations about the static Bose condensate are small, and thus we ignore the second term in (|34|), as 
discussed above. The momentum distribution of atoms in the Bose-condensed gas is then 

n c (k)= I dre lk r $5(r) / dr' e - lk r '$ (r'). (35) 



As discussed in Section 0, the order parameter of a ID Bose-condensed gas is given by the harmonic oscillator ground 
state wavefunction [see Eq. (H)] in the radial direction, and a Thomas-Fermi order parameter along the axial direction. 
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The momentum distribution then reduces to 

n c (k) = / dre l ^ r 4> p (x,y)y/n (z) 



x / dr'e- <kr V p (x',yOA/^0, (36) 



where the static condensate density no(z) is defined in (^|) Q. 
The integrals in (J3q) may be evaluated analytically to give 



n c (^,fc p y=4^^e-'?^^^M, (37) 
.91 d ^ 

where J\ is the Bessel function of the first kind of order 1. Using Eq. ( |37| ) for the momentum distribution of a Bose 
condensate, the IA dynamic structure factor in Eq. (plh for a Bose-condensed gas is given by 



7T /!() l2hu> r Ji(Zoy/m/2huj r (LJ ~ L0 r )) 

2 giD V m (lu — oj r ) 2 



Sf A (q = q t i,u>) = ~Z "V. /: 2 ^. (38) 



To determine the momentum distribution for a quasi-condensate, we must first evaluate the single-particle correla- 
tion function in Eq. (E^) in the static limit (t = 0), 



Cv t v-( r ' r '; < = °) = ^p( x ' v) ( I ) p( x 'i y')\/ n o( z ) n o( z ') 

xe -iF a( ,, z '\ (39) 
where F s (z, z') = F(z, z'; t = 0). Using the result in Eq. ( g8| ) for the phase fluctuation operator, one obtains || 

3=1 

x[P J (z/Z )-P J (z'/Z )} 2 

x [1 + 2N(huj j )] , (40) 

where N(Hu)j) is the Bose distribution function. In approximating this expression further, we make use of two facts. 
The vacuum fluctuations are small [0 and secondly, we are interested in the finite temperature regime where the 
thermally induced low energy phase fluctuations destroy the long-range order of the Bose condensate. We may then 
use the quasi-classical approximation in which the temperature is larger than the important low energy excitations. 
Replacing N(fauij) by kBT/huij, Eq. ( ff"o| ) reduces to 

f.m = 4 -^y 2j + 1 



ZT d hw z j(j + 1) 
x[P J (z/Z )-P ] (z'/Z )f. (41) 



Here the T — excitation energies are Uj = lu z \J j{j + l)/2 and following Ref . H , we have introduced the "degeneracy 
temperature" Td = Nhjj z . The j-sum over the modes may be evaluated exactly, 



F s (z,z ) = 



3 T d frw z 

(l-z/Z )(l + z'/Z ) 



In 



(l + z/Zo)(l-z'/Z ) 



(42) 



We note that Eq. (^) agrees precisely with that obtained in Ref. || by solving the coupled Bogoliubov equations, as 
it should. 

From Eq. (^) , it follows that the temperature Ts at which phase fluctuations destroy the long-range phase coherence 
of the Bose condensate can be characterized by || 

T, = (43) 
Mo 



8 



Working to lowest order in the separation \z — z'\ in Eq. (fl2|), the single-particle correlation function can be shown 
to decay as exp(— \z — z'l/l^) near the centre of the trap, where the correlation length is || 

h = Z ^. (44) 
A quasi- ID Bose gas in a highly elongated trap has a phase fluctuation temperature given by M 



32 /x 



with the correlation length 



rp3D 

l% D = Z Q -±r. (46) 



These results are essentially the same as those for trapped ID Bose gases which we work with. 

We can now use these results to calculate the momentum distribution for a quasi-condensate (denoted by the 
subscript qc), 



n qc (k) = G^(-k,k) = J dve ik r J dv'e-^ r 'C i , H {v,v'-t = 0). 



(47) 



The integration over the radial coordinates may be done analytically, and we find that the momentum distribution 
of particles in the radial direction is a Gaussian, 



n qc {k x ,k y ,k z ) = ^l 2 p e- l e {k * +k ~y ] / dze lk * z / efe'e -4 *** 

J — oo J— oo 

x^/n (z)n (zi)e~i F ^> z '\ (48) 

The dependence of n qc on the axial momentum k z is clearly complicated. Substituting the expressions for no (z) in 
Eq. (§) and F s {z,z') in Eq. Q in to Eq. (gjL the axial quasi-condensate momentum distribution must be evaluated 
numerically. Substituting Eq. (48]) into Eq. (3l]), the I A quasi-condensate dynamic structure factor is 



S Ta{^ = <lzi,u) = — \ hzr— n qc{kz = ~ ^r)), (49) 

z7r y lnw r 



where uj r is the recoil frequency and 

/oo />oo 
dze^ z / dz'e- lk ^' ^/n Q (z)n (z')e-^ F ^ z '\ (50) 
-co J — oo 

In Fig. [l], we compare Sia for a true Bose condensate and a quasi-condensate [see Eq. (ff9|)] as a function of the 
frequency u), measured relative to the recoil frequency tu r . The width of Sia for the quasi-condensate is significantly 
wider than that of a true Bose condensate, as expected. The resulting axial momentum distribution has Doppler width 
Aia proportional to l/Z^, which increases (since decreases) with increasing temperature. For a true condensate, 
the Doppler width may be evaluated analytically. Using Eq. (^8|) for the IA dynamic structure factor, the Doppler 
width using Eq. (|3^) is 

A} A = V2^-. (51) 

This Doppler width is inversely proportional to the condensate Thomas-Fermi axial length Zq, which for a true Bose 
condensate is the coherence length. Aside from the numerical factor \/2, the width is the same as given in Ref. pl| for 
a condensate in an anisotropic 3D trap. The fact that the Doppler width is inversely proportional to the condensate 
length Zq is intuitively expected. The Doppler width in Eq. (|5l]) is a measure of the spread of the momentum 
distribution of the condensate, which is inversely proportional to the spatial size of the condensate Z . 

We can also calculate the Doppler width of the IA dynamic structure factor of a quasi-condensate, as given by 
Eq. ([l9]) . We numerically fit this quasi-condensate IA dynamic structure factor to the Gaussian expansion of Eq. j32j ) 
at a series of temperatures to obtain A^. The results are given in Fig. 0, where we plot (A c IA /A q j C A ) versus T/T^. 
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Since the quantity (Aj A / Aj c ^) is proportional to (q z l^/mZo), Fig. || shows the decrease of the correlation length 
of a quasi-condensate as a function of temperature T. 

In a very recent paper, Gerbier et. al. [ [3l| have given an alternative evaluation of the static phase fluctuation 
average F s {z 1 z') using a local density approximation (LDA). This has the important advantage that one can obtain 
a simple analytic expression for F s (z,z') in Eq. (|4l|). They find that F s (z,z') is proportional to \z — z'\/l c f > (Z), but 
now the correlation length 1$ depends on the centre of mass axial position variable Z = h(z + z'), 

l^Z) = l^l - (Z/Z0) 2 } 2 (52) 

In Ref. [|T), it is shown that this simple LDA calculation of F s (z, z') is numerically an excellent approximation to the 
expression analogous to Eq. (|4l]) for a highly elongated trap for T > 8T#. This is in the temperature range of interest. 
The analytic expressions developed by Gerbier et. al. J3l| for F s (z, z') have been successfully used in analyzing their 
experimental Bragg scattering data || in terms of the altered momentum distribution of a phase-fluctuating ID 
quasi-condensate. 



IV. DYNAMIC SINGLE-PARTICLE CORRELATION FUNCTION 



We now want to extend the discussion of the well-known results given in Section III and consider the dynamic 
correlation function defined in Eq. (|29|). This is the major new result of the present paper. Using Eq. (Pq), the 
exponential argument in Eq. (^9|) is 

j=i ■> 
x[Pf(z/Z ) + P 2 (z / /Z ) 
-2P j (z/Z )P j (z'/Z ) cob(« 3 -*)] 

x[l + 2N(ftv j )]. (53) 
We can separate the static and dynamic contributions in a natural manner, 

F{z,z';t) =F s (z,z')+F d (z,z';t), (54) 
where the static contribution is given by Eq. (^2|) and the dynamic contribution is defined by 

Fd(z ' Z ' ;<) = 2 Zo{tt% S J§^Y ) P ^/Zo)P J (z'/Z ) [1 - cosfot)] . (55) 

The dynamic correlation function is then given by 

C^ H {v,r';t) = C w ( r)I ';i= 0) e -i F '^*'- t \ (56) 
which has the associated single-particle spectral density 

C v , +V) (r,r';a;) = ^t v ,(r,r / ;t = 0) / dte^e^^'^. (57) 



To evaluate the Fourier transform in Eq. (p71) , we make a "one-phonon approximation" , as discussed in Ref. 1 32 for 
a 2D harmonic crystal, and expand exp\—hF^\z, z'\ t)] to first order in F^. This approximation is good at isolating the 
low frequency (i.e., the long-time) behavior of the dynamic correlation function. In Ref. J32j, this kind of first order 
expansion was shown to give results essentially identical to the exact result in the long wavelength, low frequency 
region of interest. To lowest order, this expansion gives 

n 1 ' \ _ n 1 1 4- — n^diDksT 
l -^V t '/'v r ' r — Wit^lr, r ; t — ui — 

°° 27 + 1 

x E -rr-^Pj(z/Z )P J ^/Z ) [6(u> + u>j) + 5(uj - lu,)} . (58) 
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Since the temperature is assumed much larger than the lowest energy excitations, the energy spacing between the 
excitations is much smaller than the temperature. In this high temperature regime, we may thus approximate the 
sum by a quasi-classical integration over the excitation mode frequencies, which gives 



C^i^(v,r' ]lo) — - — ^C^^(r,T f ;t — 0) 



TTgipk B T t 
Z (tkJ z 

' dv ^ t \ Pu{z/Z )Pu{z'/Z ) [5{uj + u(y)) + 6(u - w(i/))] , (59) 
l v \ v + 1J 



where the mode frequencies = u> z v(y + l)/2 now depend on the continuous variable v. 
The integral over v in Eq. ([39|) may be evaluated using the delta function identity 

*(/(*)) = £ it^)!^-^)' ( 6 °) 

where f(x{) = and f'(xi) ^ 0. Using this in Eq. (p9|), the single-particle dynamic correlation function has the 
spectral density at low frequencies 

» — i / i \ 2\f2irgi D k B T , P v ^(z/Z )P 1/ ^(z'/Z ) 

C^{t,t ;ui) = 3 (r,r ;t = 0) ' \' 9{v(u) - 1 , 61 

where 



,( w)s -i + ^ + 2^. (62) 

The behavior of a uniform ID Bosc gas is quite different. A quasi-condensate can only form at T = 0. For non-zero 
temperatures, the system cannot be distinguished from a normal non-Bose-condensed gas. As is well known, the 
key signature of a quasi-condensate in the T = uniform ID Bose gas is that dynamic correlation functions for low 
frequencies and momentum have a power law form (see, e.g., Refs. ]l8| , p6| , |32] |) 

CWfo") ~ j^^jlZITiF . (63) 

where c = y^gno/m is the phonon velocity with no the uniform superfluid density and A" = nhy/no/mgiD is the 
Luttinger parameter. As a function of frequency w, for low momentum, the spectral density in Eq. ( |63[ ) will exhibit 
a characteristic power law singularity at uj = ck, in place of a delta function 5(u> — ck) characteristic of a 3D gas. 

The expression for the dynamic correlation function in Eq. (|6l"|) o nly describes long wavelength excitations. As 
in the case of the static correlation function discussed in Section [II, we consider the Fourier transformed dynamic 
correlation function 

C^(-k,k;w) = J dve^ r J dv'e-^'C^^vy-^). (64) 
Substituting the expression in Eq. @ into Eq. (§4|), gives 

C^(-k, k; u) = iniy^+^C^w), (65) 
where the quasi-condensate axial single-particle correlation function is defined by 

1 2\f2irg 1D k B T 



lu z Z {Hu) z y J_ 00 J_ 00 

I T~\ Tl\ -±F«(z z') P »(u)( z / Z o)Piy(ij){ z '/Z ) 

+ 1\ 



(66) 



The integrals over the axial directions z and z 1 in Eq. ( |65| ) have been evaluated numerically, with the results shown 
in Fig. | 
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For comparison, we also calculate the correlation function corresponding to an ordinary Bose condensate at T — 0. 
As in our calculations of the dynamic correlation function for a quasi-condensate, we ignore density fluctuations. The 
correlation function for a Bose condensate in a ID trap is given by [see Eqs. (§) and (§)] 

CV + </>( r > r ; *) = <M X > y)<t>p( x 'i y')^n G {z)n {z'), (67) 
so that the corresponding momentum and frequency dependent dynamic correlation function is 

G 0t ^(k,-k;w) = iirl 2 p e- l >* +k ^C c (k z ;w), (68) 
where the condensate axial single-particle correlation function is denned by 

_ 27r 3 /i JfQgqfcg) 

C (fc z ;w) = rn fl w - 69 ) 

This pure condensate result is strongly peaked about zero frequency, i.e., a delta function S(u>). 

Because of the step function in Eq. mV) , the dynamic correlation function will have no contribution from frequencies 
for which < 1, which corresponds to frequencies \u>\ less than the axial trap frequency u) z . This gap clearly arises 
because of the finite-size of the system. All frequencies below the lowest excitation energy have zero-weight in the 
dynamic correlation function, as expected. The same gap would, of course, also occur in the single-particle spectral 
density of a non-Bose condensed ID gas above Tbec- In contrast, in the case of a trapped Bose-condensed gas, there 
is still a large contribution at zero frequency arising from the condensate itself. Thus, a trapped ID quasi-condensate 
exhibits the absence of true long-range order through the presence of a gap in the correlation function spectrum 
for frequencies below that of the lowest energy collective mode (u> z ), just as a non-Bose condensed gas would. This 
effect also occurs for highly elongated traps, whose static correlation properties were considered in Ref. ||. As we 
noted earlier, aside from minor numerical differences arising from different wavefunctions for the low energy phase 
fluctuations, the behavior of a highly elongated trap is the same as the ID system we are considering in this paper. 

In the present section, we have calculated the dynamic ID single-particle correlation function directly from the 
expression in Eq. (p8|), approximating the discrete sum over the excitations by an integration. This does not give 
any simple analytic expression analogous to Eq. ( |63| ) which clearly exhibits the properties of a quasi-condensate. An 
alternative might be to use the LDA approach which Gerbier et. al. |3l| have recently used to evaluate the static 
correlation functions in a highly elongated trapped Bose gas in the high temperature limit (see discussion at the end 
of Section H|). 



OUT-COUPLING EXPERIMENT 



It has been recently shown (T^j that an out-coupling experiment, based upon two-photon Raman scattering, provides 
a way of measuring the single-particle correlation functions of a Bose-condensed gas. In this section, we apply the 
formalism of Ref. |nj (see also Refs. (ll| |l6)) to discuss how a Raman scattering experiment might be used to probe 
the dynamic single-particle correlation function of a ID quasi-condensate. So far, experimental interest has focussed 
upon the static correlation function discussed in Section III. We discuss an experiment involving out-coupling of atoms 
from a trapped hyperfine state to an untrapped hyperfine state. 

As in Ref. [ fl7| , two external laser fields produce a Raman transition from a trapped hyperfine to a second hyperfine 
state. Atoms in the trapped hypefine state are labelled by "1" and those in the untrapped atomic hyperfine state are 
labelled by "2". The energies of the trapped and untrapped hyperfine states are denoted by e\ and £2, respectively. 
The difference in energy between these two hyperfine states is denoted by Ae = £2 — 6%. The chemical potentials of 
the trapped and untrapped gases are denoted as fi\ and ^2, respectively. The difference in the chemical potentials 
is denoted as A/i = fi2 — fJ-i- Similarly, the annihilation operator "01(2) destroys an atom in the hyperfine state 1(2). 
The external potential that describes the perturbation of the laser fields can then be written as 



V ext (t) = 7 / dr 



Jq-T-iujt 



V>i(r)V4(r) +H.c. 



(70) 
u>b is the 



where 7 denotes the laser field strength, while q = k a — k(, is the wave vector difference and ui = l 
frequency difference between the two laser beams. We set h = 1 in this section. 

In Ref. fl^l , it is shown that the current of atoms flowing from a trapped gas to an untrapped gas induced by the 
external perturbation in Eq. (fOI) is given by (within linear response theory) 



/r /*OG 
dr J dr' J dte^^'le-^e^e-^-^-^ 



(r', r; -*)C^t (r', r; -t) - (r, r'; t)C^ (r, r'; t) 



(71) 
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We use the notation 

C AB (vy-t) = ^A{v,t)B(v')) (72) 

to describe the correlation function of two operators A and B in the interaction picture, with the grand-canonical 
Hamiltonian K = H — HiNi — determining the time-evolution of operators. The operator ^Vi(2) is the number 

operator for atoms in hyperfine state labelled by 1(2). 

We may rewrite Eq. ( |7l| ) f° r the out-coupling current in terms of Fourier transformed single-particle correlation 
functions [jlTj 

SI out (q, w) = 5I^ t (q, w) + «W (q, c), (73) 
where the two contributions 5I^ t and c)/,^ to the out-coupling atom current are given by 

<yi£kq>oO = 7 2 Re |dr| dr'e^(- r ') |~ ^ 

x t^t^ (r', r; -a/)C^t (r', r; J + u - Ae - A/x) 

-C^l r'^O^t^ (r, r'; -J u> + Ae + A/i)] (74) 



and 



uj" - uj' - uj + Ae + Afi 



(75) 



where the symbol P refers to principal value integration. The single-particle correlation functions are even in time 
and thus their Fourier transform is consequently real. The only component of the integrand in Eq. (|75|) for 5I^ t 



which is imaginary is the term proportional to sin[q ■ (r — r')]. Since the trapping potential is invariant under spatial 
inversions r — > — r, the correlation functions are also invariant under spatial inversions. The integrals over r and r' 
in Eq. (fq) thus vanish in the imaginary term, since sin[q • (r — r')] is odd under spatial inversion. As a result the 



contribution SI^ t is zero, and we are left with 6I$ t in Eq. Q74|). 

We may rewrite the expression for the out-coupling current in Eq. ( |74| ) in terms of the spectral density of the second 
gas A^\^ (r, r'; uj) using the well-known identity 

^('/i-w) - N 2 (uj)A^ 2 (v,v';uj), (76) 

where N 2 is the Bose distribution function for the untrapped gas evaluated at temperature T 2 . For a dilute, classical 
gas, one has N 2 (uj) <§C 1 and the expression in Eq. ([74]) reduces to 

/[• POO J f 

dr J dr'e^-^ J -^C^ (r', r; -w')^ (r, r'; J + uj - Ae - A/x). (77) 

The atoms in the second hyperfine state can be described as a uniform, non-interacting 3D gas fTij ]. Using the uniform 
gas spectral density 

A^ 2 (k,uj) = 2nS(uj-e k2 ), (78) 

where e^ 2 = k 2 /2m — [i 2 = e^ — \i 2 is the excitation energy of an out-coupled, untrapped atom, the linear response 
current SI ou t reduces to 



f dk 

5I 0Ut (cL,ui) = 7 2 Rc J — -^C^(-(k + q),(k + q); 



cj + /i -£k). (79) 



We have introduced the "detuning parameter" Co = lu — Ae and for notational convenience, we have also dropped 
the distinction between trapped and untrapped hyperfine states in Eq. (|79|). The dynamic single-particle correlation 
function and chemical potential no [see Eq. (ETi)] in Eq. ([79]) refer to the trapped gas. The result in Eq. (|79|) shows how 
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the out-coupling atomic current is formally given directly in terms of the dynamic single-particle correlation function 
of the trapped gas |l7]|. The expression in Eq. ( |79|) applies equally to a true Bose condensate and a quasi-condensate. 

We first calculate the out-coupling current arising from a ID Bose condensate at T — 0. We neglect the contribution 
from the noncondensate atoms (the thermal cloud), as in our discussion of the quasi-condensate in Section Using 
Eq. (pq) for the dynamic correlation function of a Bose condensate, the out-coupling current in Eq. (|79|) is 



« c „ t (q,^)=4^ 7 2 Re J ^ e -W k *+^ k y+^C c (k z + fo; Q + ^ £k)? (80) 

where C c is the condensate axial single-particle correlation function defined in Eq. ( |69|) . Similarly, using Eq. (|65| ) for 
the dynamic correlation function of a quasi-condensate, we find that the out-coupling current from a quasi-condensate 
is 

5I% t ( q ,Q)=4Trl^Re f ^e-^+^+^+^C^ + q z ;u + Mo - e k ), (81) 



where C qc is the quasi-condensate axial single-particle correlation function defined in Eq. (66). 

Since the radial confinement of the trapped gas is significant, the momentum distribution in the radial direction is 
very broad. The low frequency, long wavelength fluctuations in the axial direction are difficult to resolve in the total 
out-coupling current. It is better to measure the out-coupling current from atoms with low momentum. Our results 



in Section [TV] for the dynamic single-particle correlation functions in Eq. fl68| ) and Eq. (|65| ) can be used to calculate 
the momentum-resolved atom current for a quasi-condensate, and then compared to that for a true Bose condensate. 
The momentum-resolved out-coupling current from a trapped Bose condensate, using Eq. (|8C|), is 

SI c oat (q.kjC) = e^e-KiK+K) J * I *± &{ Q + Mo - £k+q ), (82) 

gi D kj 



which is strongly peaked (i.e., a delta function) about u> = ui — Ae = £k+q — Mo- hi contrast, from Eq. (|81|), the 
out-coupling current from a trapped quasi-condensate is 

W o ^(q,k;w) = 47r 7 2 Z^Ree- i ?^+ fe ?)C« c (fc 2 ; ( I) + M0-£k+ q ). (83) 

Using Eq. (|6^) for C qc , we observe that there will be no contribution to the momentum- resolved out-coupling current 
for \uj +_Uo — £k+q| < w z , in contrast to the atom current in Eq. (|82] ) from a Bose condensate. 

Fig. shows our results for the momentum-resolved out-coupling current from a true Bose condensate and a 
quasi-condensate, both broadened to take into account energy resolution. There is a gap in the quasi-condensate 
momentum-resolved current centred about £k+q ; a s there would be in the case of a quasi-lD Bose gas in a highly 
elongated cigar trap (see the discussion in Section [Tv| ). If a coherent Bose condensate is present, then there is a 
significant out-coupling current in the centre of this gap arising from the tunneling of atoms from the condensate 
(see Fig. ^). The gap in the quasi-condensate momentum-resolved out-coupling current arises simply because the 
excitations have a minimum value given by the axial trap frequency lo z . 

VI. CONCLUSIONS 

In this paper we have studied both the static and dynamic single-particle correlations in ID trapped Bose gases. 
As is well known from studies in uniform systems, the unique physics of the ID Bose gas emerges because of the 
long wavelength, low energy fluctuations of the phase, which destroy the long-range order of the Bose condensate. 
This leads to distinctive features in the single-particle correlation functions. In a uniform Bose gas, the key signature 
of a quasi-condensate is that single-particle correlations exhibit power law behavior at low frequencies and long 
wavelengths. 

Recently there has been considerable attention paid to how the static single-particle correlation functions are 
modified in ID trapped Bose gases (J, |, ^, 0, |3l[] as a result of phase fluctuations. These predictions have been 
largely confirmed by recent experiments q, 1. In the present paper, we have reviewed the calculation of the single- 
particle correlation function in ID quasi-condensates, with emphasis on the momentum distribution of atoms and 
how this can be measured using Doppler-scnsitive Bragg scattering experiments H. Our discussion of the static 



single-particle correlation functions in Section III is needed as background for the generalization in Section IV to 
time- dependent single-particle correlation functions. We also show that both the static and dynamic correlation 
functions can be evaluated using the excitations given by a quantized ID hydrodynamic model p3| , as a simpler 
alternative to previous discussions. 
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In Section M, we review a recent theory of out-coupled atom currents tunneling from trapped Bose gases [|15| |16|, [L7| . 
This kind of two-photon Raman scattering experiment might be a way of measuring the frequency dependent single- 
particle spectral density of a ID quasi-condensate. This would be a natural extension of the use of Bragg scattering 
to study the static momentum distribution of quasi-condensates || . 



APPENDIX 



The impulse approximation for the dynamic structure factor S^q, to) is often written down without much discussion. 
The physics behind the IA used in reducing Eq. (||) to Eq. (^|) is that for very large momentum and energy transfers, 
the motion of the atoms from their initial equilibrium distribution is only probed over very small distances and times. 
Thus the atom dynamics can be treated as if the atoms were free particles. Taking (h = 1) 

V>(k + q, t) = e- lBk +^(k + q, 0) (A.l) 



^(M) = e" M $ t (k.0) J 

one sees that Eq. (Eh reduces to 

—5 (w - (E k+q - E k )) (vAkMk + q)^ (k')V(k' - q; 



(A.2) 



(A.3) 



In the IA, since q is very large, one can approximate the final energy Ek+ q by the free atom kinetic energy £k+q = 
(k + q) 2 /2m. We do not necessarily know the initial energy Ek (which may be complex) of a typical atom in the 
equilibrium distribution. However, clearly one has -Ek+q 3> Ek and thus Ek is negligible. Moreover, the energy 
conservation factor in Eq. (A.3) can be reduced to 



6(u> — q 2 /2m — q • k/m), 



(A.4) 



where we only keep contributions which are consistent with the IA. Again within the IA, the static correlation function 
in Eq. (A.3) can be approximated by 



^(k)^(k + q)) (^(k')^(k' - q)) + (^(k)^(k' - q)) (i>(k + q)^(k')) • (A.5) 

The correlation functions in the first term are negligible for large q, while those in the second term may be rewritten 
to give 



S(cbw) = 



dkdk' 



5(lj — q 2 /2m — q ■ k/m) 



V' t (k)^(k'-q)) ( 5(k'-k-q) 



^(k)^(k'-q))(^(k')^(k + q) 



(A.6) 



Clearly for large q 3> k, k ', the correlation functions in the second term in Eq. (A.6) are negligible. Thus we finally 
arrive at the usual expression [B4, 33 for 5(q, u) in the IA, as given in Eq. (H) of Section O. 
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FIG. 1: The frequency dependence, relative to cj r = hq z /2m, of the dynamic structure factor in the impulse approximation 
for a Bose condensate (top panel) and a quasi-condensate (bottom panel). These results are for N = 10 4 , giD/huj z l z — 10, 
kBT^/TkUz = 7.1, lo z — 2tv x 5 Hz, lj p — 2tt x 10 4 Hz and uJ r /co z = 10 5 . 




FIG. 2: Behavior of the correlation length of a quasi-condensate (Doppler width of the dynamic structure factor) as a function 
of temperature. The parameters used are the same as in Fig. hi. 
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FIG. 4: Momentum resolved out-coupling current. The parameters TV, giD, T^, lu z and uj p are the same as in Fig. |l| with 
k x — k y — k z = 0.01/Z Z , q x = 2/7 z , q y = q z — and 7 = 0Au z . The true T = condensate current is the dashed line. The 
results have been Lorentzian broadened by T/uj z — 0.25 to reflect finite energy instrumental resolution. 



